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Mueller and Jones matrices have been thoroughly studied as mathematical tools to describe the manipulation 
of the polarization state of classical light. In particular, the most general physical transformation on the 
polarization state has been represented as an ensemble of Jones matrices, as Vi^V^ . But this has generally 
been directly assumed with out proof by most authors. In this paper, we derive this expression from simple 
physical principles and the matrix theory of positive maps. © 2013 Optical Society of America 
OCIS codes: 030.1670, 030.6600, 070.5040, 260.5430. 
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The mathematical description of the polarization state 
of a light beam dates back over a century and a half to 
the work of Stokes [1], who introduced the four parame- 
ters which bear his name to specify the beam's polariza- 
tion state [2]. The connection between these parameters 
and the stochastic theory of electrodynamics was estab- 
lished by Wiener [3] and Wolf [4,5]. 

When an optical beam passes through an optical sys- 
tem, its state of polarization will be transformed. Tra- 
ditionally, the description of these transformations has 
relied on either the Mueller matrix formulation [6] (for 
the Stokes parameters) or the Jones calculus [7] (for the 
electric field components, or more generally, for the po- 
larization matrix, also called the coherency matrix). 

It is generally assumed on physical grounds that the 
most general transformation on the state, in the polariza- 
tion matrix formalism, can be represented as an ensemble 
of Jones matrices [8], and based on this, the properties 
of the most general Mueller matrix can be derived [9,10] 
for the Stokes formalism. The proof of this assumption 
based on rigorous mathematical properties of the two 
formalisms has been lacking. Recently, Simon et al. [11] 
have addressed this problem based on properties of pos- 
itivity. 

In this paper, we address this problem differently. We 
prove that the most general transformation is indeed 
an ensemble of Jones matrices, based on Choi's work 
on positive maps [12] and simple physical assumptions 
about the state. This is similar to Choi's theorem for 
completely positive maps [13]. 

Consider a classical beam of light propagating in the 
z direction. The magnitude of the electric field in the x 
and y direction are taken to be probabilistic ensembles 
given by complex analytic signals Ei{r,t) and E2{r,t) 
respectively, where r is the position vector. The functions 
must be analytic signals, meaning when considered as 
functions of complex t, they are each analytic in the lower 
half of the complex t plane. This ensures the Fourier 
integral synthesizing the electric field converges [5, 14]. 

The polarization state of the beam of light is given by 



the 2x2 polarization matrix $(r,t), defined as follows, 

= 1 = 1,2 (1) 

where () denotes the ensemble average over different pos- 
sible realizations, and position and time dependence have 
been suppressed. 

Alternatively, the four element Stokes vector S can be 
used to represent the polarization state [1]. It is related 
to $ as follows, 
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(2) 
(3) 



where cr*^ is the identity matrix, and cr^, cr^, and are 
the three Pauli matrices cr^ , a^, and cr*' respectively. Ein- 
stein summation notation has been used, i.e. repeated in- 
dices are summed over. Lowercase Latin letters run from 
1 to 2 (corresponding to the two Cartesian components 
of the field) , while lowercase Greek letters run from to 
3. The polarization matrix or Stokes vector contain all 
the physical information about the polarization state of 
the beam [2], and are different ways of mathematically 
representing the same information. 

When the beam interacts with a linear optical ele- 
ment, generically called a filter, its polarization state is 
transformed. The most basic type of filter, which we call 
a simple filter, linearly transforms the transverse elec- 
tric field vector via the matrix T, the well-known Jones 
matrix [7], viz., 

(4) 



E- - TijEj. 



This is equivalent to the following unitary transforma- 
tion on the polarization matrix, 



(5) 



Turning to the Stokes vector, any linear transformation 
on the state must have the form 



(6) 
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where M^i, is the 4x4 Mueller matrix [6]. Equations (2), 
(3) and (5) then imply that the Mueller matrix corre- 
sponding to a simple filter can be expressed as [2,8, 10], 

1. 
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AiT (g) T*)A~^ = -A{T (g) T*)At 



(7) 
(8) 



where (g is the tensor product, and A is the matrix whose 
rows are the vectorization of the Identity and Pauli ma- 
trices, given by 



A 



1 

1 
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(9) 



Mueller matrices of this kind are called pure Mueller ma- 
trices, or Mueller- Jones matrices [10]. 

Natural questions that arise at this point arc how does 
one represent the most general type of optical filter and 
what are the relevant Jones and Mueller matrices? 

It has been a generally held axiom, motivated by phys- 
ical intuition, that the most general kind of optical fil- 
ter is an ensemble of simple filters [8-10]. Therefore, the 
transformation on the polarization matrix is then repre- 
sented by an ensemble of Jones matrices, rather than a 
single one, as follows. 



where e is the index over the elements of the ensemble, 
and Pe is the probability of realizing the e*'* ensemble 

element {J2ePe = !)• 

From eq.(7) we find that the Mueller matrix for a gen- 
eral filter is given by. 



Y,p^-Tr[a^^T,a^Tl] 



(11) 



where the M^t^ are pure Mueller matrices, each derived 
from a single Jones matrix in the ensemble. This ex- 
pression tells us that any physically admissible Mueller 
matrix is given by a convex linear combination of pure 
Mueller matrices. That is, the set of all Mueller matrices 
is the convex hull of pure Mueller matrices 

At first glance, eq.(6) seems to suggest that all ma- 
trices M^i, that map the set of physical Stokes vectors 
into itself are physical Mueller matrices. This condition 
of mapping Stokes vectors to Stokes vectors turns out to 
be a necessary but not sufficient condition for a physi- 
cal Mueller matrix [9-11]. The other necessary condition 
has to do with the positivity of the polarization matrix, 
as shown in ref. [11]. 

In the remainder of this letter, we derive eq.(lO) from 
basic mathematical properties, making use of a theo- 
rem due to Choi, regarding positive linear maps in C*- 
algebras [12]. Let !F{^) denote the operation of a gen- 
eral linear filter on the polarization matrix $. We assume 
J"(<i>) satisfies the following axioms, 



1. The map F is linear. That is, for any set of coeffi- 
cients {cj, we have J^[Y.i ^I'^i) = Y.i CiJ^{^i)- 

2. The map F is positive. That is, is a positive 
2x2 matrix for any positive 2x2 matrix $. 

3. The simplest expression for does not include 
the conjugate of the input matrix, $*. In other 
words, F is holomorphic. 

The validity of the first axiom is based on the su- 
perposition principle, assuming of course the different 
polarization matrices correspond to independent electric 
fields (i.e. the electric fields underlying the polarization 
matrices being added are mutually uncorrelated) [15]. 
The second axiom holds true since the output of is a 
physical polarization matrix, and therefore positive. 

We now justify the third axiom. From the definition of 
$ in eq.(l), we see that applying the complex conjuga- 
tion operator to $ is equivalent to applying it to the elec- 
tric fields. That is, it is equivalent to Ej — > li = I7 2. 

Reverting back to the theory of optical coherence in 
ref. [5] , we see that the electric field of a beam at a given 
point must be given by the real component of the an- 
alytic electric field, expressed as, Ej{t) = \[E''p{t) + 
ii?j*''(t)], where j = 1,2. The real and imaginary parts 
E'f\t) and Ef{t) must fo rm a Hilbert transform pair 



(10) (conjugate pair), i.e. E^p {t) = U^E^j'' ' {t)) , where % is 
the Hilbert transform. This relationship must hold for 
the complex electric field to be an analytic signal, a nec- 
essary condition for the field to be physical. Applying 
the complex conjugation operation to E.j{t) means flip- 



ping the sign of E^^'{t), and the overall field becomes 

E*{t) = ^[E^'^\t)-^Ef\t)] 

However, E*{t) does not describe an admissible com- 
plex electric field because its imaginary part is not the 
Hilbert transform of its real part (there is a sign dis- 
parity), and therefore it does not form an analytic sig- 
nal. In other words, complex conjugation of the complex 
electric field is an unphysical operation, inadmissible by 
the underlying formalism of classical linear optics. Even 
though the real part of the signal is unchanged, it is 
rendered unphysical by the imaginary part, which deter- 
mines the relative phase of the signal, e.g. in elliptical 
polarization. The conjugation operation is similar to the 
phase-conjugate mirror in non- linear optics [16]. 

We can also consider an example where we have the 
electric field Ej{t) ~ Ej{u!)e^'^*duj for some real 

Ejito). Suppose the point in question passes through the 
element that conjugates the electric field at time t = 0, 
altering it for ah future time. That is, Ej{t) — >■ E*{t) = 
Ej{-t), for t > 0. Therefore, 



j°°^Ej{uj)e"^^duj t<0 
Eji^t) t > 



(12) 



This means the conjugation operation leads to spon- 
taneous time reversal in the signal for all future time. 
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It also implies that the system somehow retained infor- 
mation about its previous state long before encounter- 
ing the hypothetical conjugator optical element. This of 
course is not physically possible for any isolated physical 
element. 

To derive the form of J", we take the first two axioms. 
Together, they state that is a positive linear map that 
maps the space of 2 x 2 matrices onto itself. In the theory 
of C* algebras [12], the most general form for such a map 
is well known. 

It is given by, 

^($) = ^F.$F;+^M^,$*W^t^ (13) 

i 3 

where <&* is the transpose of $. Recall that $ is Hermi- 
tian as can be seen from eq.(l). Therefore its transpose 
is equal to its complex conjugate ($* = <!>*). Applying 
axiom 3 to eq.(13), we see that the term involving the 
transpose must be removed altogether, i.e. Wj = Vj. 
Then we have the final form of the most general physical 
transformation upon the polarization matrix $ as, 

F{^) = Y,V,^V^. (14) 

i 

Therefore, eq.(14) describes the most general operation 
upon a polarization matrix $. It is of the form of the 
Jones matrix ensemble in eq.(lO), setting Vi = y/plTi. 

We have proven the the validity of the expression in 
eq.(lO) as the most general physical transformation on a 
polarization matrix, using some basic mathematics and 
minimal assumptions. Wc have also given physical rea- 
sons why the transpose map is inadmissible, despite its 
preservation of positivity, equating it to the unphysical 
(and acausal) conjugate map. This treatment will break 
down in higher dimensions, since eq.(13) only applies 
in the case of 2 x 2 matrices. Also note that eq.(lO) 
and eq.(14) are of the form of a completely positive 
linear map (a stronger condition than a positive linear 
map [13]). 

In the field of quantum information, quantum channels 
are transformations that have this same form for some 
unrelated reasons involving tensor products of different 
spaces [17]. 
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